UNCLASSIFIED 


AD NUMBER 
NEW LIMITATION CHANGE 
TO 


Approved for public release, distribution 
unlimited 


FROM 
Distribution authorized to U.S. Gov’t. 
agencies and their contractors; 
Administrative/Operational Use; JAN 1953. 
Other requests shall be referred to 
Defense Nuclear Energy, Washington, DC. 


AUTHORITY 


DNA/IMTI memo, 5 May 1974 


THIS PAGE IS UNCLASSIFIED 


UNCLASSIFIED 


AD NUMBER 
CLASSIFICATION CHANGES 
TO 


unclassified 


restricted 


AUTHORITY 


23 Oct 1958, per doc markings 


THIS PAGE IS UNCLASSIFIED 


ppsansanaras aia A ai aaa ac Ba a AALS RTARTA BRAM 


pnuPepereeeaaneatenen nase 


SN ee ee ee ee tee 


Wd | This document consists of 69 pages 
@ S/ A = e No.. 2 i) 8 of 357 copies, Series A 


 (Dur-soF Secpiuie. * seal 


a iat DAMAGE TO TREES— 
- ISOLATED CONIFERS. 


REPORT TO THE TEST DIRECTOR 


. a 

~~ A. A. Brown, Peeyeateitver 
; | R. K. Arnold, 

W. L. Fons) 

-F. M. Sauer 

W. E. Reifsnyder 


Ger, | Meg 
Caer cassie ; 


Classification canceiled (or changed to 


by authority f 2 7- 7o0a7 
| ot, date /O~ KAT se 


bh. 
Forest Service, U. S. Department of Agriculture C 


rite 


| | 
ff a es 


PET AEE OO ECT AVL ORL pH PEST Salt HAN te Siar nA A allel i 


<7 oo 


RESTRICTED DATA 


way 


tatiana Aha I ie tet Linn, sienna iia teak ok ok eel ote ih et hacia aha! in eniaaemh e 
a pete pet SS 


CONFIDENTIAL 
Secarity taformation 


ABSTRACT 


Project 3.3, aiming at prediction of blast damage to forests from 
atomic explosions, sought to establish the following relationships for 
isolated coniferous trees: (1) Tree damage in terms of stem breakage, 
branch breakage, and defoliation in relation to peak overpressure and 
positive phase duration time; (2) deflection-time and strain-time 
functions for tree stems following arrival of shock wave. 


Four pine trees and a weighted die. on an aiuminun ogee 
lollipop) were exposed at each of four stations, Trees approximately 
45 ft high and 12 in. in diameter were cut from the Nevada, National 
Forest and placed in concrete foundations at the test site. Lollipops, 
32 in. in diameter and weighing 380 1b, were mounted 14 ft above the 
ground on 4-in. aluminum I-beams, which in turn were’fixed in concrete 
foundations. Strain meter-oscillograph systems recorded etrain-time 
relationships on one tree and the lollipop at gach station. Motion 
pictures were used to scale deflection-time gata, (0n each tree scratch 
gages recorded maximum strain at base of wn and at 1-ft level, and 
snubber wires measured maximum stem deflection at a selected reference 


int . a 
po Ba 


a” 
_ Six trees were broken op Shot 3 and one of the remaining on Shot 
4. No defoliation or bra breakage occurred; tree crowns broke up as 
a unit with a stem bresk usually near the base of the crown. This tree 
breakage was relatec with two ratios: shock peak dynamic pressure to 
dynamic pressure for breakage due to steady wind; and tree period to 
positive phase duration time. Lollipop deflection-time history was 
also related to the latter ratic and to tue ratio of maximum deflection 
to static geflection aesociated with peak dynamic pressure. 


Good agreement was found hetween Shot 3 lollipop deflections and 
those theoretically predicted using measured overpressures. Breakage 
and deflection data on Shot 4 indicate that actual peak dynamic 
preésure was greater than that calculated from peak overpressure. This 
obgfervation agrees with Sandia Laboratory dynamic pressure measurements 
ani Los Alamos Scientific Laboratory smoke velocity data. 

\ 

Blast damage to forests from atomic explosions is primarily a 
function of aerodynamic drag associated with particie velocity. Peak 
dynamic pressure is a more compatible parameter than peak overpressure 
for predicting blast damage to trees and other structures which have 
comparatively long periods and are susceptible to damage by aerodynamic 
drag. The lower limit of peak dynamic pressure for complete breakage 
of isolated conifers is estimated to be in the order of 0.7 psi when 
associated with l-sec positive phase duration, This value of dynamic 
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4 pressure corresponds to a sea level overpressure of 5.5 psi for Mach | . 
yl reflection. Some breakage will occur at lower peak dynamic pressures tl es 
ty because of strength verteitone: in tree stems, cas 


& peri e —— oe 
} Natural forests” or ‘prepared forest stands can be instrumented a 
j economically to study blast damage from atomic explosions with no re- v } 
‘ quirement for outside power or timing signals. 
It ie recommended that more static aaa as be made to es- 

tablish variations in static strength of growing tree stems, and that 

future effects tests consist of a graded series of shots to increase 

the range of duration times as well as overpressures. 
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CHAPTER 1 


OBIECTIVE 


Project 3.3 was one phase in a research program aimed at pre- 
diction of blast damage to forests from atomic explosions, Studies 
connected with Project 3.3 were confined to isolated coniferous trees 
in order to simplify the problem and to establish the following re- 
lationships, which provide a field check on blast damage prediction: 


1. Tree damage in terms of stem breakage, branch breakage, and 
defoliation in relation to peak overpressure and positive phase dura- 
tion time, 


2. Deflection-time and strain-time functions for tree stems fol- 
lowing arrival of shock wave. 


To provide a further -heck on computational systems used for 
damage prediction, deflection- and strain-time relations were measured 
on a simple spring-mass sy-tem which was designed to represent an 
"ideal-reproducible" tree, 


These studies of isolated coniferovs stems were preliminary to 
more comprehensive projects which deal with tree stands of varying 
densities, shapes, and combinations of tree species. Therefcre it was 
important for Project 3.3 to develop and field-test methods and tech- 
niques of measuring motion and strain on tree stems subject to blast 
from atomic explosions in order to instrument a forest stand at reason- 
able cost. 


9 
CONFIDENTIAL RESTRICTiO DATA 
fecerty etrartes 


ATOMIC ENERGY ACT 1945 


when, 
Sn ee 


CONFIDENTIAL 


Sacerity tatermation 


CHAPTER 2 


BACKGROUND AND THEORETICAL DATA 


2.1 GENERAL 


Preliminary analysis)/ of damage trees by shock wind indicated 
that at sea Level, conifers 24-in. DB with 50 per cent crown and a 
natural period of 2.5 sec would be expected to break eat ap overpressure 
of approximately 7.5 psi in the Mach stem region. The purpose of this 
preliminary work was more to develop methods or analysis and delineate 
problem parameters than to establish a radius of tree breakage. This 
analysis provided the basis for several field studies which followed. 


Studies of aerodynamic drag in tree crowns2/ indicated the nature 
of forces which blast winds from atomic explosions might exert on iso- 
lated tree crowns. Variation in tree crown drag was found to be due 
primarily to bending upon application of drag force. Dry weight of 
crown appeared to be the best parameter with which draz could be corre- 
lated. 


Bending and breakage characteristics of ponderosa pine stems sub- 
jected to static loading have been analyzed ,4/ and indicate that maxi- 
mum atress at the point of static failure is of the order of 3,000 
lb/sq in. Trees subjected to static loading near the center of 
pressure of the crown generally broke near the base of crown, 


VY Operations Research Office. Preliminary Study of the Consequences 
of an Atomic Explosion Qver a Forest. ORO-T-10%. Washington, 1950. 
102 pp. 

2/ Diemeter at breast height, 4-1/2 rt above the ground, outside bark. 
3/ U. S. Dept. of Agriculture, Forest Service, Division of Fire 
Research. Experimental Investigation of Aerodynazic Drag in Tree 
Crowns Exposed to Steady Wind—Conifers. Phase xeport for Cperaticns 
Research Office. December 20, 1951. 19 pp. 

4/ U. S. Dept. of agriculture, Forest Service, Division of Fire 
Research. [ree Breakage Characteristics Under Static Loading--Ponder- 
98a, . Report in preparation for Armed Forces Special Weepons 
Prdject. 
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2.2 COMPUTATION SYSTEM. 


Further analysis of tree--stem breakage by shock wind, 2/ wnich 
employed data from work described above, established significant 
parameters of tree characteristics and shock wave characteristics which 
affect the probability of breakage and bic:-down following atomic ex- 
plosions. 


In this analysis the tree was replaced by a single-mass, spring- 
mass system with the mass placed at the position of the center of 
pressure of the crown drag force. Natural period of the simplified 
system ecualed that of the tree, and the mass had the drag associated 
with the crown. Restoring force characteristics of the spring were 
determined from static deflection measurements made on the tree stem. 
The viscous damping constant was determined from damped free vibration 
of the stem. 


The equation cf motion equates drag, restoring, and daming forces 
to the inertial force 


a as = px Dy) - R(y) - dy (2.1) 


where V is the relative velocity between the shock particle velocity 
and the velocity of the mass center, y is the displacement of the mass 
center, and ,O is the air density. Density and particle velocity are 
assumed to have a linear function with time, “quation 2.1 wag then 
numerically integrated in the manner outlined by Timoshenko. Figures 
2.1 and 2.2, results of calculations mate to dave, were usad to locate 
stations for Project 3.3. Similar calculations were made on v2ighted 
beams to predict maximum deflections for each shot. These calculations 
were used to determine strength requirements and strain recorder gain 
settings. A more detailed vresentatian of the method is in Appendix A. 


5) U. S. Dept. of Agriculture, Forest Service, Division of Fre 
Research. Analysis of Tree-Stem Breakage by Shock wind--Ponderosa 


Pine. Report in preparation for armed Forces Special weapons Pro ject. 


6/ Timoshenko, S. Vibration Problems jn Engineering. (2nd ed.), New 
York: D. YanNostrand Co., Inc., 1937. p. 126. 
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Fig. 2.1 Calculated Breakage Overpressures for Isolated Ponderosa Pine. 
Based on Sea Level Free Air Overpressures for a 20 KT Bomb at 
2,000 Ft Burst Height. 
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Fig. 2.2 Normalized Breakage Shock Pressure Ratio for Isolated Ponder- 
osa Pine, Static Shock Pressure Ratio Based on Shock Peak 
Dynamic Pressure Equal to That Associated With Breakage Due 
to Steady Wind as Calculated From Normal Shuck Relationships. 
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CHAPTER 3 


PREPARATIONS AND METHODS 


3.1 STATIONS AND TREES 


Four trees and a "lollipop")/ were placed at each of four stations 
on the Forest Service line in Area 7. Stations covered a range of 1.9 
to 5.3 psi peak overpressure for Shots 3 and 4 as shown in Table 4.1, 
Figures 3.1, 3.2, 3.3, and 3.4 are station photographs. The plot plan 
and station layout are shown in Fig. 3.5. 


[Or ren cote emer ee mre ctmnee conn n a g men RSET TN gp Rena On Sore copies = memaaRpAy To sare erm “4 


. . 
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Fig. 3.1 Station I Showing Lollipop and Trees A, 38, C, and D (Left 
to Right) With Crowns 40, 60, 80, and 30 Per Cent of 
Total Height. View From Direction of Ground Zero. 


Trees, approximately 45 ft high and 12 in, DBH, were cut from 

Site Class 4 and 5 areas on the Charleston Ranger District of the 
Nevada National Forest. Whole trees were hauled by special wrailer to 

. the Nevada Proving Grounds and anchored by placing butt-ends 6 ft into 
prepared concrete foundations, Stems were then grouted in place and 
branches trimmed to form typical classes of tree crowns, which can be 
seen in station photographs. Actual cutting was done 20 days before 
Shot 3. Trees were positioned at Stations I and II before Shct 2 and 
at Stations III and IV 10 days before Shor 3. 


1/ Weighted disk on an aluminum beam. 
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of Total Height. View From Direction of Ground Zero. 
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Fig. 3.3 Station III Showing Lollipop and Trees A, B, C, and D 
(Left to Right) with Crowns 50, 50, 50, and 50 Per Cent ( 
of Total Height. View From Direction of Ground Zero, | 
ee % 4 | 
, @ ae ‘| { 
a rn | een 7 sik aaron. ai. aye aire ood 
e Fig. 3.4 Station IV Showing Lollipop and Trees A, 5, C, and D 
(Left to Right) With Crowns 50, 60, 80, and 4C Per Cent 
of Total Height. View From Direction of Ground Zero. 
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of the crown (which approximates center of pressure), and natural peri- 
ods of vibration were measured prior to Shot 3. 


| “ Reference points were established on each tree stem at the center 
Following each shot physical characterisiice of broken trees were 

analyzed by measuring stem diameter inside bark at 5-ft intervals and 
by determining weight and moisture content of foliage and branchwood. 

: After Shot 4, natural periods of standing trees were taken again. 

i Static load-deflection and load-strain relations were then determined 

5 for "A" trees by pulling each tree with a cable attached to the stem 

at the reference point and measuring deflection, strain, and load until 

the tree broke, Physical characteristics of these and the remaining 

standing trees were then analyzed. These measurements and relations. 

are summarized in Tables 3.1 and 3.2 and in Fig. 3.6. 


TABLE 3.1 


Mechanical Characteristics of Trees 


Deflection Spring Constant Strain Constant 
for Lb/Ft In./In,./Ft 
Breakage*™ 


Reference | Center of Reference | Center of 
Po Pressure Point Pressure 


int 
146 237 627 672 
289 512 835 1250 
150 218 710 
183 2% 938 


* Deflection for breakage as a result of static loading applied at 
reference point. 
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TABLE 3.2 
" Phyuical Characteristice of Trees 


Dry Weight 
of Foiiage 


Nominal 

Per Cent or 
Tree Height 
-. and Branches 
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® All trees were western yellow pine, Pinus ponderosa var scopylorum 
Engelmann, except III D and IV C which were red fir, Abies magnifica 


Murray, 

b Diameter at breast height, 4-1/2 ft above ground, outside bark. 

© Diameter inside bark 

@ Height at which deflection measurements were made end static loads 
applied. 
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Fig. 3.6 Static Strain, Load, Deflection Relationships 
for Instrumented Trees 
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3-2 LOLLIPOPS 


A singie-mass system was designed as an “ideal-reproducible tree" 
with a natural period typical of trees exposed. Natural period and 
other characteristics of the single mass system (Table 3.3) were re- 

7 produced in lollipops, which were constructed by mounting a 32-in. 
diameter steel, concrete-filled disk on an aluminum I-bean. 


These I-beams were anchored in concrete foundations in the same 
manner as the trees and were guyed perpendicular to the radius from 


ground zero. 
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Fig. 3.7 Lollipop at Station I as Seen From Direction of Ground 
' Zero, Steel Guard at Base Protects Strain Gage. 


| 2/ 2.72-1b, 4-in, I-beaw, alloy 61 ST 6. 
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“TABLE 3.3 
Lollipop Characteristics 


Spring Constant - 341 1b/f% 


Natural Period 1.18 sec 


Disk Weight - 380 lbs Disk Area - 804 sq in, 


Beam Weight - 34 lbs Beam Area ~ 351 eq in, 


Equivalent Total Rquivalent Total 
Disk Weight - 387 lbs Disk Area - 1005 aq in. 


0.75 


Strain Constant - 1920 #in/i:,* Base Fixity? 
tt 


Height to Center of Disk ~ 14 fe 


. prin. /in. « strain in microinches per inch 
b Ratio of actual to theoretical spring constant values 


3.3 INSTRUMENTATION 


3.3.1 Strain-time Measurement 


Baldwin SR-4 strain gages were mounted 7 in. above the base 
on the front and rear of each lollipop. A "strain meter" (Fig. 3.8) 
was mounted on the ground zero (tension) side of tree "A" at each sta- 
tion, 1 ft above the base, so that the gage could be affixed to the 
outermost wood surface (Fig. 3.9). Gage assemblies were protected 
from blast damage by foam rubber padding and steel guards similar to 
that shown in Fig. 3.7. 


Strain meter and strain gage bridge output were recorded on 
Sanborn Model 127 oscillograph recorders equipped with Model 140 strain 
gage amplifiers, The system was independent of external power supply 
or timing signals, A block diagram of the power and relay system is 
shown in Fig. 3.10, 


Power was supplied by four 1l2-volt, 50 amp. hr automotive 
wet batteries connected in series-pavallel. Mach oscillograph power 
circuit contained a separate 24- to 110-voit AC converter, Power con- 
sumption during the shot was only 4 amp. hrs so that battery voltage 
drop was small, 
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Fig. 3.9 Strein Meter Mounted an Tree 
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Fig. 3.10 Block Diagram — Cscillograph Power and Relay Systems 
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All relays, except the photocell relay, were 24-volt aircraft 
mechanical-electrical latch in type as insurance agsinst maifunction 
due to shock. 


The two alarm clocks were set at H-10 hrs so that alarms would 
ring at H~1 hr and actuate the parallel toggle switches to close the 
Agastat circuit. Strain amplifiers were thus powered during a one-hour 
warmup perioc for stabiiization. The fireball tripped the photocell 
(Fig. 3.11), actuating the photocell relay, which in turn closed the 
two tape drive motor relays and cpened the Agastat operational relay. 
After a time delay of 2 min., all power was turned off by the Agastat. 


> + 5 7 8 % 
ae : .O ~N oN Oe . < 
— Re ae. cade ow, race AA, pw eS ae res i & E ara) e es oe Pars oo 


Fig. 3.11 Photocell and Mount Above Instrument Shelter 


Upon recovery the toggles were opened and all relays were reset 
with an eiectrical reset switch. Alarms were reset and the system was 
immediavely ready for operation. 


The pnotcell cireuitd/ was too sensitive, This defect was cor- 
rected by spraying the tubes with aluminum paint so that a Wabash Press 
No. 50 flash bulb would trip the relay at 2 ft but not at 10 ft. 


All circuits leading to strain meters, strain cages, amd photo- 
cells were shielded microphone cable, In addition, the first tvo 
cables ran through conduit providing magnetic as well as electrical 
shielding, All microphone cables and oscillographs were grounded to 
the conduit. 


ee ee oe 


3/ Thomas 1. Monahan et al. The Effect of Thermal Radiation on Materi- 
als. Naval Material Laboratory, New York Naval Shipyard. A.F.S.@.P, : 
Report WT-311, Operation BUSTER, p. 26. i 
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3.3.2 Maximum Strain Measurement 


Baldwin deForest scratch gages were mounted to measure maxi- 
mum strain at the following locations on trees: 


1 f% above ground away from ground zero - all trees mf 
1 ft above ground toward ground zerc - trees B, C, D 
Base of crown away from ground zero - all trees | 


3.3.3 Maximum Deflection Measuycnent 


Deflection stubbers were made by attaching piano wire to 
each tree at the reference point and to the center of the lollipop disk 
on the side towards ground zero. These wires led on a 45° angle to 
emall snubbing blocks staked down on the ground zero side of the tree. 
Snubber wires were pulled tight and the length of long tails extending 
beyond the snubbers was measured prior to each shot. Snubber-wire ends if 
were then remeasure¢ upon recovery. 


3.3.4 Motion Pictures 


Camera stations for motion pictures were designed to place 
the plane of the film parallel to the radius from ground zero in order 
to measure deflection versus time, Figure 5.12 shows the camera view of 
tree movement, Also shown are aluminum foil bands which were placed on 
each tree 10 ft and 15 ft above ground level for scaling purposes. A 
foil band was also placed at the reference point. Thirty-five milli- 
meter motion pictures were taken by the Army Pictorial Service Division, 
Office of the Chief Signal Officer, TUMBLER-SNAPPER Project 9.1. 


L. a | : Re aan, a 
~ a eat 


Fig. 3.12 Station I Viewed From Motion Picture Camera Station 
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CHAPTER 4 


RESULIS 


4.1 GENERAL 
Overpressures and positive phase duration times to which trees and 
lollipops were exposed are shown in Table 4.1. 
TABLE 4.1 


Overpressures and Positive Phase Duration Times at 
Project 3.3 Stations, TUMBLER Shots 2, 3, and 4* 


Overpressure (PSI) Positive Phase 
10-Ft Level Duration Tines (Sec) 
Station Peak Incident Ground Level 
Shot Shot 2 
I 2.1 0.42 0.98 0.98 
II 1.6 0.43 1.02 1,02 
IV = = 1,09 1.09 


* Interpolated from data furnished by the Armed Forces Special Weapons 
Projecte 9 May 1952 


Six trees were broken on Shot 3, and one of the remaining trees on 
Shot 4 ss shown in Table 4.2, No main branches were broken by either 
blast. There were a few scattered twigs on the ground at each station 
after Shot 3 but none after Shot 4. 


Unbroken trees and lollipops did not appear to be damaged struc- 
turally. Strain and deflection calibrations and natural period measure- 
ments made before Shot 3 and after Shot 4 agreed within limits of 
experimental error, 
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TABLE 4.2 
a 
Characteristics and Circumstances of Tree Breakage 
Station Over- | Height of 
and Shot | pressure* | Break Above Type of Break 
(PSI) Base (Ft) 
3 5.3 1 Diagonal bending failure 
3 5.3 14 Diagonal bending failure origi- 
nating at knot on tension 
3 4.3 2 Short bending failure 
3 4.3 9 Brash failure at knot cluster 
3 4.3 19 Brash failure at one knot on 
tension 
3 3.6 3 Diagonal tension failure 
{ 
4 4.1 7 Diagonal tension failure 
* Ten-ft-Level values, se Table 4.1. " 
ey mia q 
Cee ' 
ee : | 
L. < “ . i Sa ren me ath . 4 
‘ the. _ am : 
pert Ts 
@ ° Pee . a cae 
a 
3 oe | 
Fig. 4.1 Breakage of Trees B, C, D at Station II, Shot. 3 
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Fig. 4.2 Typical Stem Failures. A - Brash Failure at 
One Kiiot on Tension; Tree II D. B - Diagonal 
Tension Failure; Tree iII D. 


Tree breakage and non-breakage data for Shots 3 and 4, Figs. 4.3 
and 4.4, were normalized in the manner described in Appendix A. Dy- 
namic pressure ratio was used rather than shock-pressure ratio because 
the former is compatible with both regular and Mach reflection types of 
shock when loading is principally due to drag. Dynamic pressures for 
breakage due to steady wind were estimated from tree characteristics 
given in Table 3.2, Conditional breakage regions in Figs. 4.3 and 4.4 
correspond to that predicted theoretically (Fig. 2.2) when it was moved 
parallel to itself so that it separated non-breakage and breakage data. 
No weight was given to the three trees which failed at knots and which 
consequently fell in the region of non~breakage in Fig. 4.3. 


4.2 DEFLECTION, TIME HISTORY —~ TREES 


Deflection-time histories and maximum deflections were scaled 
from motion pictures whenever clarity permitted, Since a comparison 
of motion picture time of shock arrival with actual shock arrival times 
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showed wide variation from nomin.1 film speed of 64 frames/sec, film 
speed was estimated from osciilograph records of instrumented trees, 
Table 4.3. Oscillograph tape seed timed before and after each shot 
showed a masimum disparity of 8 per cent. 


TABLE 4.3 


Oscillograph Time of Arrival and Film Speed Data 


Oscillograph | Film Speed (Frame/Sec) 
b ‘Based On 
Shot | Station Time of Oscillograp 
Shock Arrival] Tape Speed° 
; 


2 4l. 
53.4 = 
3 57,2 54.7 
: 57.5 a 
84.5 60.8 
; 29,1 31.8 
62.4 64.6 
4 64.3 64.6 
16.4 73.5 


. Interpolated from data furnished by the Armed Forces Special Weapons 
Projecte 9 May 1952 


b Not to be taken as true times of arrival -- for comparison only. 
© During time of tree motion. 


| : Deflection-time histories and maximum deflectians were determined 
from strain meter data (Fig. 4.5) by dividing measured strain at 1 ft 
by the static strain constant for the reference poirt, 
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Fig. 4.5 Tracing of Strain Meter Record, Tree I A, Shot 3. 
Strain = 1090 # in./in,/division; Time = 2.68 


division/sec. 
Deflection-time histories computed by the method described in 
Section 2.2 are compared with those scaled from motion pictures and 


those determined from strain meter data in Figs. 4.6, 4.7, and 4.8 
Maximum deflection data are compared in Table 4.4. 


sree pMETER 


00 0.2 G4 Ow 
TIME PROM SHOCK ARRIVAL (SEC) 


Fig. 4.6 Motion of Reference Point Under Blast 
Loading, Tree I A, Shot 2 i 
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DEFLECTION (FT) 


Qo 
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TIME FROM SHOCK ARRIVAL ( SEC) 


Fig. 4.7 Motion of Reference Point Under Blast Loading—Trees, Shot 3. 
See Table 3.2 for Reference Point Height. 
- Motion Picture Data 


—-—~-~ Strain Meter Data 
———— Calculated 
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Fig. 4.8 Motion of Reference Point Under Blast Loading—Trees, Shot 4. 
See Table 3.2 for Reference Point Height. 
&§ ———--— Motion Picture Data 
a Strain Meter Data 
s—-~-— Calculated 2 
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a Questionable--snubber wire moved, 


TABLE 4.4 
Maximum Tree Deflectiona 


& Deflection for breakage 


at 

3 
A 

2.4 


© Measured strain at 1 ft divided by statio strain constant for reference point. 
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4.3 MAXIMIM DEFLECTION —~ LuiLIPOPS 


Maximum deflections were measured from motion picture films, but 
relative movements of lollipop centers of pressure were too small for 
weasurement of deflection-time histories. Oscillograph data were pro- 
cessed in the same manmer as tree strain meter data. Maximum deflec- 
tions for centers of pressure are shown in Table 4.5 and are compared ; ot 
on ea dimensionless basis in Fig. 4.9. 


TABLE 4.5 


Maximum Lollipop Deflections 


Deflection 
(Sec) 


Sa 


BERS BIS8a 


J 


COPN OF PE 


* 7 — data not reliable due to twisting of lollipop v8 


4.4 SCRATCH GAGE STRAIN MEASUREMENT 


Maximum strain values in Table 4.6 were read with a toolmaker's 
microscope to the nearest u.0002 in. (equivalent to a strain of 0.1 
fin./in.) directly from scratch gage targets. igure 4.10 presents a 
typical target record. Only 6 out of 70 gages failed to record. . 


Data t'rrom simultaneous free vibration records taxen on Tree IV A 
with scratch gage and strain meter on opposite sides of the stem are 
compared in Fig. 4.11. 
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Fig. 4.10 Scratch Gage Record and Diagrammatic Explanation 
(Magnification = X 32), Tree IV B, Shot 3, 1 Ft 
Above Ground Away From Ground Zero. 
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STRAIN METER READING (10344 INCH/INCH) 
Fig. 4.11 Comparison of Strain teasurement — Scratch Gage 
Versus Strain Meter 
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TABLE 4.6 


Meximum Strain Compariso. — Scratch Gages and Strain 
Meters, TUMBLER-SNAPPER Shots 3 and 4 


a 


= 
Maximum Strain, 10? 4 in,/in,| Maximum Strain, 10 in,/in, 
1 | Rc | 


IA 7,55° 21 21.4 17.3 
IB 8.4 32.5 19.3 33.3 
1¢ 5.2 11.6 M M 
ID U 1.9 M HM 
II A 6,49 10.3 7.9 9.1 
II B Be 10.9 M M 
1c B 12.9 M M 
II D U 12.7 M M 
III A 5,08 9.8 5.7 8.0 
III B 5.0 8.4 6.5 7.1 
: III ¢ U 5.5 4.9 44 
III D 7.6 15 M 


IVA 2035 6.2 4.4 5.7 
IV B 5.6 8.0 3.8 3.6 
Ivc 5.7 5.9 3.3 2.9 

6.9 3.7 4.0 


. Gage location on tree: 


F-1, tension, 1 ft above ground facing ground sero 
R-1, compression, 1 ft above ground facing eway from ground sero 
R-C, compression, at base of crown facing away from ground zero 
> Underlined values recurded by strain meter located in F-l position 
© U = Record unreadable or missing 
d M «= No record; tree broken by Shot 3 
@ 5B = No record; tree break near gags 


FO AER tate eR R SO AML R MER 


senetee 


a 


37 
CONFIDENTIAL RESTRICTED DATA 
Security inderastion 


ne Sie eteeatanneemantementntbeieiod 


, 
cae en es ae Tes a cet ge a 
a 


5.1 GENERAL 


Predictions of minimum overpressures which will break forest trees 
ere based on calculations which assume that stem loading is the result 
of aerodynamic drag of the tree crown due to after-shock winds, Aero- 
dynamic drag characteristics of crowns have been determined experi~ 
mentally and correlated on a dimensionless basis with physical crown 
characteristics. Calculations to date have been based on linear time 
variation of after-shock particle velocity and density or on the theo- 
retical time history appropriate to free air as shown in Fig. A.2. 
However, the validity of applying free air time histories in regions of 
regular and Mach reflection remains to be verified. 


$.1.1 plat £ Shot ta 


Maximum particle velocity and density following shock — 
peak dynamic pressure — determines the magnitude of forces acting on 
trees. In the Mach reflection region peak dynamic pressure was calcu- 
lated from peak overpressure by normal shock equations; in the region 
of regular reflection the calculation was based on incident and re- 
flected shock overpressures by using oblique shock equations. Result- 
ant peak dynamic pressure values are shown in Table 4.5 and correlated 
lollipop deflection data for Shot 3 are shown in Fig. 4.9. 


If the fact that Stations I and II were in the regular 
reflecticn region on Shot 31/ (Stations III and IV were in the Mach 
reflection region) is neglected and peak dynamic pressure at those two 
atations is computed by normal shock equations, this failure to account 
for oblicue shock produces a 35 per cent discrepancy between correlated 
lollipop deflections at the first two stations and the last two stations, 


These correlation studies indicate that peak dynamic pres- 
sure and duration time are the fundamental shock parameters for pre- 
dicting or comparing damages to structures such as trees, suspension 
bridges and towers, which have comparatively long periods and are af- 
fected primarily by aerodynamic drag. These parameters are particularly 
applicable for damage comparison between regular and Mach reflection 
regions, 


V Test Data Rurnished by the Armed Forces Special Weapone Project. 
15 May 1952. 
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5.1.2 Comparison of Shot 3 and Shot 4 Results 
Differences between results of Shot 3 and Shot 4 are sun- 
marized below: 
, 1. Lollipop maximum deflections at Station I were higher ( 


on Shot 4 than on Shot 3 even though peak dynamic 
pressure determined from reported overpressures for 
Shot 4 was less than for Shot 3, and the positive phase 
duration times were the same (Tables 4.1 and 4.5). 


2. ILellipop maximum deflections for all stations on Shot 
4 were higher than predicted when this prediction was 
based on results of Shot 3 and overpressures as re- 
ported for Shot 4. 


3. Maximum deflection from snubber data and maximum base 
strain from strain meter and scratch gage data for 
tree I A, Shot 4, were higher than for Shot 3. 


4. Tree I B, unbroken on Shot 3, was broken on Shot 4, 
This tree appeared structurally and physicaliy un- 
changed after Shot 3. 


These differences indicate that results from Shots 3 and 4 could not be 
correlated satisfactorily even though allowance was made for Mach and 
regular reflection occurrence as discussed in the previous section. 

The only area of agreement between Shots 3 and 4 lies in the time to 
maximum deflection for the lollipops. Little doubt remains, however, : 
that higher comparative deflections of Shot 4 over Shot 3 were a true : 
phenomenon, These data tend to support the hypothesis that factors ; 
heretofore not considered caused peak cynamic pressures on Shot 4 to 

be higher than indicated by reported peak overpressure#. 


Sandia Laboratory pitot tube results2/ lend supporting 

evidence to the above hypothesis. Measured dynamic pressures on Shot 
4 at the 13-ft level were 20 to 30 per cent higher than those computed 
from measured overpressures in the 3 to 8 psi overpressure region, 

- Measured Shot 3 dynamic pressures were in substantial agreement with 
computed values, Further, but inconclusive, support ori tes in Los 
Alamos Scientific Laboratory Shot 4 smoke velocity data. Measure- 
ments at the 300-ft level indicate peak particle velocities 10 to 15 
per cent higher than those calculated from overpressure measurements 
made nearer the ground. 


NCA Sige PUBMED Poe Wtfas ont, 


2/ Sandia Laboratory. Shock~Gguge Evaluation Tests, A.F.S.W.P. Oper- 
ation TUMBLER~SHAPPER Report WT-505, Projects i9.1lce and 19.1d. 
Verbal communication with F. B. Porzel of Los Alamos scientific 
poratory, September 17, 1952. 
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| “4 Both deta are subject to scatter due to low resolution at low 

ri - overpressure level; however, they agree in order of magnitude and in- 
| 4 dicate the possibility of « 30 per cent increase in dynamic pressure 
over computed values. A 30 per cent correction added to dynamic pres- 
i sures used for correlation of Shot 4 lollipcp deflection and tree 
breakage data would bring them into juxtaposition with Shot 3 data. 

i 


An additional pressure pulse occurred on Shot 4 during the posi- 
tive phase.4/ Since this occurred after maximum digplacements of 
lollipops and trees it may be presumed that this secondary shock had 
no effect on these displacements, 


3.2 TREE BREAKAGE 


Due to the limited number of tree installations and to the fact 
that no two trees are exactly alike, prediction of tree breakage must 
lean heavily on calculation procedures and extrapolation of data by 
theoretical means. In addition, because of the reduction of KT on Shot 
3 from the originally intended bomb size, none of the fully instrumen- 
ted trees broke. These situations precluded full exploitation of the 
theory of breakage ieveloped in Chapter 2; however, breakage, non- 
breakage, and conditional breakage regions of Figs. 4.3 and 4.4 have { 
been defined. Results of Shots 3 and 4 are displayed separately since < | 
the discussion in Section 5.1.2 indicated shock phenomena were differ- 
ent in each case. 


The conditional breakage region determined from theoretical con- 
siderations alone (Fig. 2.2) is subject to certain inherent inaccura- 
cles in magnitude, Although the equation of motion may be integrated 
with sufficient accuracy, dynamic breakage deflections are in doubt 
due to a lack of complete knowledge of the increase in strength of 
green timber under high loading rates, Further, more recent static 
breakage data indicate that approximately 50 per cent of all trees will 
be stronger than the trees used as tne basis of the conditional break~ 
age region displayed in Fig. 2.2. The procedure followed in the 
theoretical analysis should produce the correct variation with respect 
to period ratio; however, dynamic pressure magnitudes must be ad justed 
in accord with field experiments as was done in Figs. 4.3 and 4.4. 


4.4 because of variations in tree characteristics and because it is 
unlikely that a tree is located just at the breaking point on any one 


ry, U. S. Naval Ordnance Laboratory. Free Air and Ground Level Pres- 
gure Measurements. A.F.S.W.P. Operation TUMBLER Report, Projects 1.3 ~ 
1.5. 15 September 1952. 
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shot. However, it is believed that the three breakage points on Shot 3 
which fall within the region of non-breakage (Trees I D, II C, II D) 
are not chance variations since these and only these trees broke at 
knots as shown in Table 4.2. It is possible that stress concentrations 
occur at certain knot combinations, More actual breakage data for 
static and dyncmic loading are needed to esteblish the validity of this 
hypothesis. 


In Figs. 4.3 and 4.4 conditional breakage regions were located to 
separate breakage and ncn-breakage date except for the case described 
above. The scatter _n these diagrams as well as the Shot 3 variations 
indicate need for more test data which, added to that already obtained, 
can be used to establish breakage probabilities associated with varia- 
tions in tree and shock parameters. 


Section 5.1.2 and the above discussion indicate that the conditim- 
al breakage region determined from Shot 3 data (Fig. 4.3) represents 
the better basis for the lower limit of dynamic pressure approaching 
' gomplete breakage of isolated conifers, 


The combination of tree periods and positive phase times of Shots 
3 and 4 grouped data at the lower end of the period ratio scale (Figs, 
4.3 and 4.4). Since periods of growing trees of the type used for 
Operation SNAPPER range from 2 to 4 sec and periods of some species 
range up to 14 sec, future teats should cover a greater range of period 
ratio, With improved methods of installation it may be possible to in- 
crease tree periods to 2 sec; however, an increase in range of positive 
phase duration times as well as overpressure should also be provided 
through variations in bomb size, 


Stem breakage occurred near the dase of the crown or below, and 
in no case did the crown break up. However, since these trees were 
selected because they had no defect and no dead limbs, breaks at such 
defective points in svems and crowns can be anticipated in a natural 
forest. 


9.3 DEFLECTION-TIME HISTORY — TREES 


Deflection-time data determined from strain measurements generally 
agreed with those determined from motion pictures when strain constant 
at small deflection was used, However, had the full static strain- 
deflection curve (Fig. 3.6) been used, excessive deflection values 
would have been obtained from strain meter data, and in some cases such 
deflection determinations would have been in oxcess of breakage, 
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It has been esteblished’/ that under conditions of high loading | 
rate, the yield point for wood increases above its static value. Suf- 
ficient data are not available in the literature to predict this in- ; i 
crease for growing trees. Project 3.3 results indicate that this | A 
increase for green tiinber is much greater thar would be predicted from | 
dry wood data. One example is found in the .act that Tree II A failed - Hl 
under static lcading at a deflection of 6.6 ft, while motion pictures 
shoved its deflection on Shot 3 to be 8.8 ft without breaking. This 4 
phenomenon cannot be explained by moisture content difference after 
Shot 3 (90 per cent) and efter Shot 4 (60 per cent), because within 
this range experimental evidence indicates that there is no appreciable 
change in the static strength of wood. 


A detailed examination of corresponding motion picture deflections 
and strain meter readings should lead to the determination of the 
stress-strain relationship under dynamic loading. Unfortunately this 
comparison cannot be made in the region close to breakage, since none 
of the instrumented trees failed. Furthermore, film speed of the 
motion pictures was in doubt, Consequently, film speed was arbitrarily 
made to agree with oscillograph tape speed by matching times at initial 
and first zero defiection positions. Agreement or disagreement between 
times to maximum deflection is coincidental. 


As a result of the high yield point when loading rate is high, de- 
flection calculations on instrumented trees were made using the linear 
load-deflection relationship as determined from the spring constant : 
instead of using the static load-deflection relationship described in 
Section 3.1. Had the latter method been used, greater deflections would 
have been calculated, and in some cases breakage would have been pre~ 
dicted. 


Drag determinations were made from previous worko/ which estab- 
lished weight ratio of foliage to branch wood and crown weight as the 
most important trce parameters affecting drag. Since SNAPPER test 
trees were cut prior to the current year's needle grovth, the weight 
ratio based on these previous studies, which were made in late summer 
when needle growth was complete, was reduced by one-third. 


Calculated maxiuwm deflections were within 20 per cent of measured 
values, Agreement among calculated deflection-time histories is con- 
sidered good in view of the number of variables involved. This agree- 
ment substantiates use of the single mass approximation used in previous 
conifer analyses and use of the free air time dependency of particle 
velocity and density as a suitable approximation in the regular reflec- 
tion and Mach stem regions. 


5/ L.J. Marquardt and T.R.C. Wilson. Strength and Related Properties 


of Woods Grown in the United States, U. S. Dept. of Agriculture Tech. 
Wasnington: U.S, Govt. Printing Office, 1935. p. 60. 


. 479, 
g. "Forest Service, gerouynamic Drag in Tree Crowns, op, cit. 
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5.4 DERLECIION-TIME HISTORY = LOLLIPOPS 


_ Derivation of the equation of motion outlined in Appendix A indi- 
cates that normalized maximm deflection is a fumction of four dimen- 
Sionless variables for any particular particle velucity and density 
history. 


3B. F_,.%s ,%, Ph, 
y, ~ ‘lor’ qr’ 1 2) 


Prior to TUMBLER~SNAPPER a series of numerical integrations was 
perforwed using a linear time dependency for particle velocity and 
density, These calculatiors covered the range of variables expected 
for the lollipops during the test as shown in Table 5.1. 


TABLE 5.1 


Range of Variavles Covered by Deflection Calculations 


0.00175 


0.10 


In the range indiceted «bove only the period ratio was significant, 
shock overpressure was second order, and the remaining two variables 
insignificant. TUMBLER lollipop data remained within the above range 
of calculations. 


A second set of integrations was perfoimed using free air decay 
of icle veiocity and density shown in Ng. A.2 for mean valuea of 
Ty/U]T, and Yg/L and pressure ratio range experienced at TUMBLER. 


Figure 4.9 compares results of these two computations, Use co? 
linear decay gives somewhat higher deflection ratior at higher period 
ratio, but the disparity is not great. 
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Maximum deflection and time to reach maximum deflection data were 
normalised as described above. Points plotted in Fig. 4.9 are based on 
averages of the three methods of measuring deflection except at closer 
stations, where twisting of lollipop beams invalidated strain gage 
records. In these cases motion picture and smubber data were used to 
calculate plotted points, 


Shot 4 data fall below the predicted curve for reasons discussed 
previously. Since the strain record went off scale on the one Shot 2 
record, estimated maximum strain and the resultant plotted point may be 
unreliable, Shot 3 results are considered to be the most reliable, and 
agreement between these experimental points and theory is good. 


Substantiation of the calculation method by these lollipop data 
and the close approximation of the linear velocity decay calculations 
to the free air. calculations indicate that previous tree breakage 
calculations snown in Fig, 2.2 which were based on e linear decay are 
not considered to be seriously in error. The trend indicated by these 
calculations should be reasonably correct; liowever, magnitudes must be 
ad freted by experiment as was done in Figs. 4.3 and 4.4. 


5.5 SCRATCH GAGE DATA’ 


Seratch gage installations were made to provide basic data on 
maxinnm svrain at base of tree stem and at base of crown to provide 
estimates of breakage-strain under impulsive loading, and to study the 
applicability of low-cost, strain~measuring instruments. 


Preliminary analysis of scratch gage data indicates that: (1) 
Large deflections resulted in considerable distortion of the strain 
curve across the l-ft section as would be expected from curved beam 
theory, and (2) greater strain occurred at the base of crown which adds 
to the evidence that breakage will occur near this point, 


Figure 4.11 establishes one-to-one correspondence between strain 
meter and scratch gage readings, These data were taken at low de- 
flection from a record similar to Fig. 4.10 and are maximum tension to 
maximum compression reacings. The maximum strain level zepresented is 
therefore 3 4 in./in., and the readings are not influenced by non- 
linearity of strain across the section of the tree stem, which occurs 
at high deflections. 


Accuracy of scratch gage readings is dependent primarily on the 
establishment of zero reading, i.e., zero strain level prisr to shock 
arrival. As shown in Fig. 4.10 this line is rather wide. Aithough 
readings were made to the nearest 0.1 4 in,/in, of strain, estimated 
accuracy on the basis of Fig. 4.11 difficulty in est.blishing zero 
is of the order of 0.5 4 in,/in,. 
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5.6 INSTRUMENTATION 
; 5.6.1 Saubbers 
Snubbers are the easiest and least costly way to measure 
. maximum deflection. Field data can be read quickly and accurately, 


Reduction of these data to maximum deflection, however, requires that 

the arc traveled by the reference point be known, This difficulty has 

been overcome in recent static breakage tests by attaching two snubbers 

at the reference point, and thus determining maximum deflection direct- 

ly from field measurements. Snubbers can be used to instrument an 

entire stand of trees, natural or artificial, which will be exposed to 
j atomic explosions. 


5.6.2 Scratch Gages 


Scratch gages represent an expedient methoa of recording 
maximum strain whenever time sequence is not important. Accuracy is 
approximately 1/10 that of strain meters but in view of low cost and 
small installation time ($4.00 and 1/2 hr per installation), this can 
ve tolerated when a large amount of data is desired. Reduction of 
data requircs a microscope and is rather tedious. Readable scratches 

: are not always obtained (only 6 out of 70 used by Project 3.3 failed to 
record), but experience in setting ecratch gage arm pressure reduces 
record loss to a minimum, 


5.6.3 Qscillograph-Strain System 


Oscillograph recording of strain meter and strain gage cir- 
cuits provides accurate measurement of time variations in strain and 
deflection, This system was well-adapted to field use and suffered no 
damage from induced signals, shock, or sand, Amplifier gain settings 
vere stable, Attachment of strain meters to trees was satisfectory as 
evidenced by the fact that records returned to their original pre-shock 
sero readings. Strain-time records may be obtained from any point in I 
the tree stem by use of this system suitably protected against blast : 
damage. } 


The principal deterrent to large-scale instrumentation to ; 
obtain strain-time history on many trees is cost. A complete oscillo- 
] graph-strain system to record strain-time at one point may be dupli- 

cated at a cost of approximately $900 and requires one day to install 
after the shelter and conduits are in place. 


5.6.4 Motion Pictures 


Motion pictures provided general visual data on tree and 
lollipop movement following shock arrival but were only partially 
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satisfactory for searurement of deflection. Disparity between nominal 
film speed and actual film speed made determination of time sequence 
44fficnlt. Hed film speed been constant and known, the disparity con- 
cerning the dynamic pressure on Shot 4 could have been resolved from 
measurenents of the displacement of vapor originating from the tree 
foliage. Dust obscured all motion at Station I and motion after the 
first maximum at other stations on Shot 3. Collapse of a temporary 
shelter cover at Station III on Shot 3 covered the camera lens before 
maximum deflections occurred. All prints of Shot 3 were fuzzy, which 
mace it difficult to fullow precisely the position of tree stems, 


5.6.3 Lollipops 


Lollipops represeut instruments which respond to peak 
dynamic pressure and may prove valuable in determining this quantity 
within forest stands, Lollipops at Operation SNAPPER failed to respond 
as intended because of their twisting motion during large deflections, 
This difficulty can be resolved by selection of a beam cross-section of 
higher torsional radius of gyra‘ion or by a redesign on tne principle 
of the ballistic pendulum, Consideration should be given to the possi- 
bility of making the lollipop portable to permit changes in its loca- 
tion between shots of a greded series. 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


6.1 CONCLUSIONS 


1. Blast damage to forests from atomic explosions is primarily a 
function of aerodynamic drag of the particle velocity. This basic 
assumption was verified by results of Project 3.3. 


2. A method of predicting tree breakage has been established and 
in part substantiated. This method relates tree breakage with ¢vwo 
ratios: shock peak dynamic pressure to dynamic pressure for breskage 
due to steady wind; and the ratio of tree period to positive phase 
duration time. 


3. Peak dynamic pressure is a more compatible parameter than peak 
overpressure for predicting atomic explosion blast damage to trees and 
structures which have comparatively long periods and are suscey “ible 
to damage by aerodynamic drag. Peak dynamic pressure criteria -ppear 
to be compatible with breakage and deflection results obtained in Mach 
and regular reflection regions, whereas peak overpressure alone does 
not permit correlation of results. 


4. Coniferous tree crowns break up as a unit with stem breakage 
usually occurring near the base of the crown. Strain on the tree stem 


is greater at the base of crown than at the 1-ft level. 


5. Peculiarities in breakage and deflection data on Shot 4 in- 
dicate that peak dynamic prc .oures experienced were greater than those 
calculated from measured peak overpressures, 


6. Agreement between measured and calculated lollipop deflectious 
for Shot 3 substantiate calculation methods used to predict tree motion 
and the use of free air particle velocity and density time history «es 
approximations for Mach end regular reflection regions. 


7. Applicability of a single-mass system to analytically predict 
tree motion following shock was confirmed by Project 3.3 results. 


8. The lower limit of shock peak dynamic pressure for complete 
breakage of ‘solated conifers is estimated to be in the order of 0.7 


psi when associated with l-sec positive phase curation. This value of 
dynamic pressure corresponds to a sea level overpressure of 5.5 pai for 
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Mach reflection. Some breakage will occur at lower dynamic pressures 

because of statistical veriations of static strength found in ap- : 
parently sound tree steas and alao because of stem Jefect such as rot 

or fire scar. 


9. Lollipops represent en instrument capable of measuring the 
total impulse due to dynemic pressure of chock particle velocity. 


10. Yetural forests or prepared forest stands can be instrumented 
economically *o study blest damage from atomic explosions with no re- 
quis went for outside power or timing signals. "he retrain zeter systen 
acturrtoly records time history of strain, snuovers provide adequate 
meesurenent of ma~imum deflection, and scratch gages measure maxima 
strain, 


6.2 BECQMENDATIONS 
6.2.1 Piture Atom Bomb Field Tests 


i. The period ratio (tree period to duration time) covered 
in Project 3.3 should be extended by using trees with longer periods 
ani by increasing the range of duration time through a graded series of 
shots. 


2. ‘'neident and reflected overpressures for the reguler 
reflection reg should continue to be reported in the 2- to 8&psi 
overpressure region. Values should be measured approximately 20 ft 
avove the ground. 


3. Dynamic pressure should be measured in preference to 
overvressure whenever such data are used for correlatior against 
ctructurai cemage due primarily to aerodynamic dreg. 


4. Investigation of shock wave disturbance due tc low 
burst height should be continued, particularly in relation to peek 
par jcle ve.ccity and density. 


*, Instrumentation: 


(a) Portable lollipops should be designed for use in 
and around forest stands. 


(b) Seratch gages should be used tu investigate 
maximum strain along tree stens 


6. Before a prepared stand is exposed to atomic explosions 
the proble= of what minimum area will constitute a semi-infinite stand 


mast be solved. 
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6.2.2 Other Work 


* 1. More static breakage tests are needed to establish vari- 
ation in static strength of growing tree stems. 


7 2. laboratory tests should be made to determine deflection- 
strain relations to the point of failure for tree stems. 
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APPENDIX A a 
GALCULATION METHODS 


A.l RQUATION OF MOTION = LOLLIPOPS 


Lollipop motion under blast loading is assumed to be that of a 
simple spring mass system acting under external force due to aero~ 
dynemic drag of the after-shock particle velocity. Internal damping 
force ie assumed negligible in compar:on to other forces, The system 
is shown schematically in Fig. A.1l. 


Fig. A.l Equivalent Mechanical System — Lollipop 
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From Newton's principle the equation of moticoa is 
my + kp (y)= (2x ODA). (a.1) 
where particle velocity relative to center of pressure velocity is given 
by 
Ve U4 - 2 y coed | (A.2) 


ant (2E CDA) is the drag component in the direction of motion, For 


ewall deflections the center of pressure follows a circular path and to 
the first approximation the beam assumes a sine curve, The angles are 
then 


62X ~My 
, om aa 5 : (a.3) 


Referring to Fig. A.1 and combining Equetions A.l and A.2 


2 2 -.2 ‘ 
ay + wig 7 SL om (B) {() ‘® oa cos 5} 


D onstaieeaiatet slices ooeniieieraamiantciaaetniieatie een tte em. A ENTE OPA ETD RARER NC 
% 


(a.4) 
coe (2+%) cos (X-2') 
Letting 
2 
=Z1,r%= t s = B 
eee eg) ae 
_ | 
2 ; 
re uty | 
Equation A.4 can be written in dimensionless form: | 
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ae +8) cos( a- 2) 


y UT, (A.5) 


For free air © «= f(x) and 2-0 = fo(x) as shown in Fig. A.2. 
U} Pi-Po 


The density ratio zay be found from the equation 


21 - LEP {1 - £5(x) | A.6 
-’ 1 { anh rs 


Equations A.5 were integrated numerically for various values of 


TF Y Yg , and to establish maximum values z, Yp , and the 
WT F) hit , ’ ’ Ye ) 


time to reach this value, Ty , for a linear particle velocity and 
density Ty, 


and for free ai particle velocity and density relationships. The re- 
sults of these calctations are shown in Fig. 4.9. 
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PARTICLE VELOCITY 
MAXIMUM PARTICLE VELOCITY 


U 
— 
t 
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%, POSITWE PHASE CURATION TIME 


Fig. A.2 Theoretical Particle Velocity and Density Variation With Time 
for Free Air, 2.2 < (Py - Po) < 8.4 PSI. Equations Represent 
Bepirical Fit Through Points Computed From L.A.S.L. IBM Prob- 
lem M Calculations Supplied by Lt. Col. F. B. Porsel, Los 
Alamos Scientific Laboratory. Subscript o Refers to Preshock 
Values, Subdecript 1 to Peak Values. 
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A.2.1 Regular Reflection 


For incident-reflected shock configurations (Fig. A.3) the 
density ratio may be found from 


Lo iw (A.7) 
22 
Po l+ # 1+ # gt 
V 


eince thermodynamic state equations are invariant under translation, 


DIRECTION OF SHOCK MOTION <== 


q, 
REGION | 


Fig. A.3 Incident-Reflected Shock in Stationary Frame of Reference 


1/ R. Courant and K. 0. Friedrichs. oubepsgnis Ege and Shock Waves. 
nterscience Publishers, Inc., New York, 1945. p. ae 
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aa 


be brought to rest by im- 
(Fig. A.3) where 


(4.8) 


Angle of incidenc ence, f8 » wae determined from burst height and station 


ground distance. 


equation’/ 


F(s) 0-74 -£@) 0) 


where 
‘. 2 
(=) ~ A, 
fy Py fl 
and 
—> ee —- 
qq = U1, + 


The peak particle velocity may then be found from the energy 


(a.9) 


In order to make the calculation method entirely clear, a 


numerical example of dynamic pressure determination is presented. 


sider the illustrative set of data: 


burst height 
ground distance 


Ibid., PP. 148-149. 


= 2500 ft 
= 3500 ft 
Gg = 1100 ft/sec 


3/ This procedure vas only approximete for Station II, Shot 3, where 


Mach refiection had formed, but the 
level. For that pressure level /3 - 


this approximation was reasonably correct. 


4/ Courant and Friedrichs, op,cit., p. 302. 
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riple point was below the 10-ft 
extr = 0.6°, which indicates that 


RESTRETED DATA 
ATOMS ENERGY ACT 1048 
L 
ry eaten il 
eS 
"2 lage 
- : 5 ws . & 
Lae ea ish deat aN 


Con- 


: fon ot tpn nen, snot + nate necerage cnee ar ; 
F ; : ose <a a. nena 
yess sie Rt I nicl A: oat } is 4 Se a ae ctr nee ery emer te 
oset Si a . pov Te aalilg ute Rae te : : prigietet.s ei ee 
pihs.. ee Had a te Spe Peat + ae ee re eee acta PT . ie 


Aa tan) 2 q 54,50 
pf - : for Ys 1.40 


From Equation (A.7) 


1 
fiw 222338 15,0 . 1.26 
Po 6 +55%9 6 + 15-0 


From Equation (A.8) 


23,0 
I 2 1 1+6 ie 0 
(2) " (e14)2 7 mat 


From Equation (A.9) 
ay 1.709 + (1 - Hd x rhe) = 2.215 
Sy 13.07 1.2645 °° % 


Hence 


U, * (u = 99\ cy = (1.104 - 1.307)1100 = -224 ft/sec 
fo 8% 


The dynami: pressure becomes 
2 2 
asl. & age (a). FA w(2) 


= 0.70 x 1.26 x 13.0 (0.203) = 0.473 pat 
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A.2.2 Mach Reflection 


” wl a normal shock, Equations A.7 and A.9 can be written i 
explicit : 


~- 


+p | 
om Se (a.10) 
Sy (+ jo) (ye? n FL) . 


A.3 — 


To reduce experimental deflection data to dimensionless form, base eat 
fixity was computed from the spring constant based on the lollipop peri- i 
od, as given by Equation A.12, and the theoretical spring constant com 

puted from the beam constants, ass a fixed end cantilever and 

loading at the center of the disk. 


2 : 
ty = A (wate + 25 Yoean) (4.12) 


F = r constant from 2 (A.13) 
computed spring constant for fixed end cantilever : 


z actual d etl 
t theoretical deflection (A.14) | 


The theoretical spring constant k under steady wind loading was 


computed ysing a distributed drag loading along the beam oy equal 
to Linn? cone plus the drag due to the lollipop dick, 2}1 ¢ 


6/ S. Timoshenko. Vibration Problems in Mngineering. (2nd ed.) , 
D, Van Nostrand Company, Inc., New York, 1937. p. 86. r 
57 | 
CONFIDENTIAL RESTRICTED DATA : 
Secarity intwmsten ATOMIC ENERGY ACT 1808 | 
; OOS gg sai 
. Sige ~. 


CONFIDENTIAL 


Sarcurity haderexstion 


Drag coefficient was assumed to correspond to the steady state valuel/ 
of 1.1 for the disk and 1.8 for the beam. The equivalent disk area 


then equals 


1¢ canter of nressure deflection due to beam drag ) 
deq * aa ( center of pressure deflection due to disk drag 


(A.15) 
Actual static deflection corrected for base fixity becomes 
012 
¥.< opr Cyheg (a.16) 


A.4 _ 


Trees are assumed to be simple spring mass systems with the same 
characteristics ag previously described for lollipops except that the 
drag term is based on weight of dry crown and is a function of displace- 
nent .8/ An additional assumption is that the crown drag is independent 
of the angle of attack. Internal damping is neglected as experimental- 
ly determined values of the damping coefficient are negligible when 
compared with the crown drag force, 


Calculations were carried out in a dimensional manner using Eque- 
tion A.4 with the peak particle velocity, density, and positive phase 
duration time appropriate for each shot and station. 


A.5 REDUCTION OF BREAKAGE DATA — TREES 


Dimensionless breakage relationships presented in Figs. 4.3 and 
4.4 may be inferred from the following analysis. From results of lolli- 
pop calculations the maximum deflection of a simple spring mass system 
may be represented to the first order as 


%. (2 
Ta Ty 


7/ Lionel S. Marks. Mechanical Engineers’ Handbook, McGrev-Hil! Book 


Company, New York, 1941. p. 1543. 
8/ Forest Service, Analysis of Tree-Stem Breakage by Shock Wind— 


Ponderosa Pine, op, cit. 
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Consider a completely linear system for which breakage occurs at a peak 
dynamic pressure of /2)U ): A static deflection may be defines on 
2 


this basis: 
(ay = AX) 


Under steady wind loading, breakage occurs at a dynamic pressure de- 
fined by 


% - (AE), 


If under impulsive loading breakage occurs at the same deflection es 
umer static loeding by steady wind, then Y, mist equal Yp. Therefore, 


(Ys)p = Yt (Z) 
Hence 
ey 
( ) is ) (4.17) 


oe) MR 
*), 


For tree breakage the system is not linear to breakage. Neverthe- 
lees, the dimensionless relationship inferred in Equation A.17 may be 
applied intuitively provided it is realized that the functional re- 
lationship may not be the same, 


VW 
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| 


be 


ES 
; 
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NOMENCLATURE 


Ag = lollipop disk avea, sq ft 

Ae = aquivalent lollipop disk area, sq ft 

Ag = lollipop beam area per ft of length, sq ft/ft 

Cg = sound speed at prevailing atmospheric conditions, ft/sec 
C,; = sound speed immediately behind shock wave, ft/sec 

Cp = drag coefficient, dimensionless 


drag function, sq ft 


i=) 
~ 
<< 
~~ 
“ 


e = 2,718..., base of natural logarithms, dimensionless 


<-- 1n  tli AAA TION ARTI ann a 


f = a functional relation 
F = base fixity, dimensionless 
g = gravitational constant, ft/sec” 
k = epring constant for ho-iz-ntal loading, lb/ft | 
ky = spring constant determined from period of vibration, 1b/ft 3 
L = height of center of mass above ground, ft 

m = mass, ib/sec?/ft 

No = flow velocity normal to incident shock, ft/sec 

Po = ambient air pressure, 1d/sq in 

Py = air pressure immediately behind shock, 1b/sq in 

P2 = air pressure immediately behind incident shock, lb/sq in 
Io = flow velocity ahead of shock, Fig. A.3, ft/sec 

4, = flow velocity behind shock, Fig. A.3, ft/sec 


restoring force function, 1b 
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t = time, sec 
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time to maximm deflection, sec 
positive phase duration, sec 
particle velocity, ft/sec 
maximum particle velocity behind shock, ft/sec 
relative velocity, ft/sec 
weight, lb 
time parameter, t/T,, dimensionless 
mass displacement along arc, ft 
fat, ft/sec 
d2y/at2, ft/sec2 
deflection of center of mass ut breakage, ft 
maximum deflection of center of mass, ft 


static deflection of center of mass associat d with peak 
dynamic pressure, ft 


deflection parameter, y/Ys, dimensionless 


dz/dx, dimensionless 
d°z/ax*, dimensionless 
an angle, Fig. A.l, radians 


angle of incident shock, radians 


an angle, Fig. A.1l, radians 

an angle, Fig. A.1, radians 

— , where 7 is the ratio of specific heat of air, 
dimensionless 

3.1416..., dimensionless 


air density, 1b sec*/rt* 
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